To appear in J. Fluid Mech. 



1 



Controlling the onset of turbulence by 

streamwise traveling waves. 
Part 2. Direct numerical simulations 

By Binh K. Lieu, Rashad Moarref and Mihailo R. Jovanovic 

Department of Electrical and Computer Engineering, University of Minnesota, Minneapolis, 

MN 55455, USA 

(Received 30 October 2009 and in revised form 30 April 2010) 

This work builds on and confirms the theoretical findings of Part 1 of this paper, Moarref) 



& Jovanovic (2010). We use direct numerical simulations of the Navier-Stokes equations 
to assess the efficacy of blowing and suction in the form of streamwise traveling waves 
for controlling the onset of turbulence in a channel flow. We highlight the effects of 
the modified base flow on the dynamics of velocity fluctuations and net power balance. 
Our simulations verify the theoretical predictions of Part 1 that the upstream traveling 
waves promote turbulence even when the uncontrolled flow stays laminar. On the other 
hand, the downstream traveling waves with parameters selected in Part 1 are capable of 
reducing the fluctuations' kinetic energy, thereby maintaining the laminar flow. In flows 
driven by a fixed pressure gradient, a positive net efficiency as large as 25 % relative to 
the uncontrolled turbulent flow can be achieved with downstream waves. Furthermore, 
we show that these waves can also relaminarize fully developed turbulent flows at low 
Reynolds numbers. We conclude that the theory developed in Part 1 for the linearized 
flow equations with uncertainty has considerable ability to predict full-scale phenomena. 



1. Introduction 

The problem of controlling channel flows using strategies that do not require measure- 
ment of the flow quantities and disturbances has recently received significant attention. 
Examples of these sensorless approaches to flow control include wall geometry deforma- 
tion such as riblets, transverse wall oscillations, and control of conductive fluids using the 
Lorentz force, to name only a few. Min et al. (2006) used direct numerical simulations 
(DNS) to show that surface blowing and suction in the form of an upstream traveling 
wave (UTW) leads to a sustained sub-laminar drag in a fully developed turbulent channel 
flow. This motivated Marusic, Joseph & Mahesh (2007) to derive a criterion for achieving 
sub-laminar drag and to compare laminar and turbulent channel flows with and without 
control. Furthermore, Hoepffner & Fukagata (2009) characterized the mechanism behind 
UTWs as a pumping rather than as a drag reduction; this is because the UTWs increase 
flux relative to the uncontrolled flow. Finally, Bewley (2009) and "Fukagata, Sugiyama & 



Kasagi ( 2009 ) independently established that for any blowing and suction boundary ac- 



tuation, the power exerted at the walls is always larger than the power saved by reducing 
drag to sub-laminar levels. This lead the authors of these two papers to conclude that 
the optimal control solution is to relaminarize the flow. 

Heretofore, sensorless flow control strategies have been designed by combining physical 
intuition with extensive numerical and experimental studies. For example, a number 
of simulations on turbulent drag reduction by means of spanwise wall oscillation was 
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conducted by Quadrio & Ricco ( 2004 ) where 37 cases of different control parameters were 
considered. Compared to the turbulent uncontrolled flow, a maximum drag reduction of 
44.7% was reported. However, analysis of the power spent by the movement of the walls 
shows that a maximum net power gain of only 7.3 % can be achieved. Even though 
DNS and experiments offer valuable insight into sensor less strategies, their utility can be 
significantly enhanced by developing a model-based framework for sensorless flow control 
design. 

This paper builds directly on the theoretical findings of Part 1 , [Moarref fc Jovanovic 



( |2010 ), where receptivity analysis was used to show that the downstream traveling waves 
(DTWs) are capable of reducing energy amplification of velocity fluctuations in a tran- 
sitional channel flow. The effectiveness of DTWs and UTWs in preventing or enhancing 
transition is examined in this work. In contrast to the current practice, we do not use 
DNS as a design tool; rather, we utilize them as a means for verification and validation 
of theoretical predictions offered in Part 1 of this study. Namely, we use DNS to confirm 
that the DTWs with parameters selected in Part 1 can control the onset of turbulence 
and achieve positive net efficiency relative to the uncontrolled flow that becomes turbu- 
lent. On the contrary, the UTWs enhance transient growth and induce turbulence even 
when the uncontrolled flow stays laminar. In spite of promoting turbulence, the UTWs 
with large amplitudes can provide sub-laminar drag coefficient. However, we show that 
this comes at the expense of poor net power balance in flows driven by a fixed pressure 
gradient. This is in agreement with |Hoepffner & Fukagata ( 2009| ), where it was shown 
that it costs more to achieve the same amount of pumping using wall-transpiration than 
pressure gradient type of actuation. Our numerical simulations show the predictive power 
of the theoretical framework developed in Part 1 and suggest that the linearized Navier- 
Stokes (NS) equations with uncertainty represent an effective control-oriented model for 
maintaining the laminar flow. 

Our presentation is organized as follows: in § [2j we present the governing equations, 
describe the numerical method used in our simulations, and outline the influence of trav- 
eling waves on control net efficiency. The evolution of three dimensional (3D) fluctuations 
around base flows induced by surface blowing and suction is studied in § |3] We further 
emphasize how velocity fluctuations affect skin-friction drag coefficient and net power 
balance. The energy amplification mechanisms are discussed in § |3.4[ where we show 
that the DTWs improve transient behavior relative to the uncontrolled flow by reducing 
the production of kinetic energy. In addition, the effect of traveling waves on coherent flow 



structures in transitional flows is visualized in § 3.5, where it is shown that the DTWs 
control the onset of turbulence by weakening the intensity of the streamwise streaks. 
In § [4j we show that the downstream waves designed in Part 1 can also relaminarize fully 
developed turbulent flows at low Reynolds numbers. We summarize our presentation and 
give an outlook for future research directions in § |5] 



2. Problem formulation and numerical method 

2.1. Governing equations 

We consider a 3D incompressible flow of a viscous Newtonian fluid in a straight channel; 
see figure [l] for geometry. The spatial coordinates (x, z) are scaled with the channel 
half height, (5, and they denote the streamwise, wall- normal, and spanwise directions, 
respectively; the velocities are scaled with the centerline velocity of the laminar parabolic 
profile, Uc; the pressure is scaled with pU^^ where p denotes the fluid density; and the 
time is scaled with the convective time scale, S/Uc. The flow is driven by a streamwise 
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downstream 

Figure 1. A channel flow with blowing and suction along the walls. 



pressure gradient and it satisfies the non-dimensional NS and continuity equations 



(u • V) u - VP + (l/i?c)Au, = V-u. 



(2.1) 



Here, Rc denotes the Reynolds number, Rc = UcS/u^ v is the kinematic viscosity, u is the 
velocity vector, P is the pressure, V is the gradient, and A is the Laplacian, A = V • V. 

In addition to the constant pressure gradient, = —2/Rc^ let the flow be subject to 
a zero-net-mass-fiux surface blowing and suction in the form of a streamwise traveling 
wave (Min et al. 2006). The base velocity, = ([/, V, ly = 0), represents the steady-state 



solution to (2.1) in the presence of the following boundary conditions 



V{y = ±1) = T2acos(cj^(x - ct)), U{±1) = Vy{±l) = W{±1) 



0, 



(2.2) 



where uj^^ c, and a denote frequency, speed, and amplitude of the streamwise traveling 
wave. Positive values of c identify a DTW, while negative values of c identify a UTW. 
In the presence of velocity fluctuations, u represents the sum of base velocity, U5, and 
velocity fluctuations, v = (u^v^w), where and w denote the fluctuations in the 

streamwise, wall-normal, and spanwise directions, respectively. 

2.2. Numerical method 

The streamwise traveling waves, considered theoretically in Part 1, are tested in DNS of 
a 3D transitional Poiseuille flow in this work. All DNS calculations are obtained using the 
code developed by'G ibsoii| ( |2007 ). A multistep semi-implicit Adams-Bashforth/Backward- 
Differentiation (AB/BDE) scheme described in Peyret (2002) is used for time discretiza- 
tion. The AB/BDE treats the linear terms implicitly and the nonlinear terms explicitly. 
A spectral method (ICanuto et alMdSS) is used for the spatial derivatives with Chebyshev 
polynomial expansion in the wall-normal direction and Fourier series expansion in the 
streamwise and spanwise directions. Aliasing errors from the evaluation of the nonlinear 
terms are removed by the 3/2-rule when the horizontal FFTs are computed. We modified 



the code to account for the streamwise traveling wave boundary conditions (2.2). 



The NS equations are integrated in time with the objective of computing fluctuations' 
kinetic energy, skin- friction drag coefficient, and net power balance, §[3j The velocity field 
is first initialized with the laminar parabolic profile in the absence of 3D fluctuations, 
§ |2.3[ this yields the 2D base flow which is induced by the fixed pressure gradient, 
Px = —2/Rc^ and the boundary conditions (2.2). In simulations of the full 3D flows (cf. 
§[3|, an initial 3D perturbation is superimposed to the base velocity, u^,. As the initial 
perturbation, we consider a random velocity field developed by Gibson (2007) which 
has the ability to trigger turbulence by exciting all the relevant Fourier and Chebyshev 
modes. This divergence- free initial condition is composed of random spectral coefficients 
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Case 


Symbol 


c 




a 






Ny 









X 








27T 


47v/3 


65 


50 


50 


1 


□ 


5 


2 


0.035 


27V 


47v/3 


65 


50 


50 


2 


o 


5 


2 


0.050 


27V 


47v/3 


65 


50 


50 


3 





5 


2 


0.125 


27V 


47v/3 


65 


50 


50 


4 


< 


-2 


0.5 


0.015 


87V 


47V 13 


65 


200 


50 


5 


V 


-2 


0.5 


0.050 


87V 


47V 13 


65 


200 


50 


6 


A 


-2 


0.5 


0.125 


87V 


47V 13 


65 


200 


50 



Table 1 . The computational domain and spatial discretization considered in simulations of the 
uncontrolled flow, DTWs with (c = 5, a;^, = 2, a = {0.035,0.050,0.125}), and UTWs with 
(c = —2, Ux = 0.5, a = {0.015, 0.050, 0.125}). Symbols identify the corresponding flow in figures 
that follow. The box sizes in the streamwise and spanwise directions are denoted by Lx and L^, 
respectively. The number of grid points in the streamwise, wall-normal, and spanwise directions 
are represented by A^i, i = {x,y,z}, respectively. 



that decay exponentially and satisfy homogenous Dirichlet boundary conditions at the 
walls. The flux and energy of the velocity fluctuations are computed at each time step. 

A fixed pressure gradient is enforced in all simulations which are initiated at Rc = 
2000; this value corresponds to the Reynolds number R^ = 63.25 based on the friction 
velocity, Ur. Owing to the fixed pressure gradient, the steady-state value of Rr is the 
same for all simulations, Rr = 63.25. In addition, we consider a streamwise box length, 
Lx = 4tt /ujx^ for all controlled flow simulations. This box length captures the streamwise 
modes kx = {0, ±co'a;/2, icj^;, ±3co'a;/2, . . .}; relative to Part 1, these modes correspond 
to the union of the fundamental {kx = {0, ^ojx^ :^2uJx-, • • •}) and subharmonic {kx = 
{±uJx/2^ ±3uJx/2^ . . .}) modes. In addition to the uncontrolled flow, we consider three 
DTWs with (c = 5, cj^ = 2, a = {0.035,0.050,0.125}), and three UTWs with (c = -2, 
ujx = 0.5, a = {0.015,0.050,0.125}). The complete list of the parameters along with the 
computational domain sizes and the number of spatial grid points is shown in table [l] The 
total integration time is ttot = 1000 5 /Uc. We have verified our simulations by making 
sure that the changes in results are negligible by increasing the number of wall-normal 
grid points to Ny = 97. 



2.3. Base flow and nominal net efficiency 
Base velocity, U5 = {U (x, t), V(x, y^ t), 0), is computed using DNS of 2D Poiseuille flow 



with Rr = 63.25 in the presence of streamwise traveling wave boundary control (2.2). 
Figured shows the mean velocity profiles, U{y) (with overline denoting the average over 
horizontal directions), in uncontrolled flow and in flows subject to selected DTWs and 
UTWs; these results agree with the results obtained using Newton's method in Part 1. 
The nominal bulk flux, which quantifies the area under U{y)^ 

Ub,n = I tu{y)dy, 



27-1 

and the nominal skin-friction drag coefficient for three UTWs and three DTWs are 
reported in table [2] For fixed pressure gradient, Px = —2/Rc^ the nominal skin- friction 
drag coefficient is inversely proportional to square of the nominal flux, i.e., 

Cf,N = -2PjUlj^. (2.3) 




(a) (b) 



Figure 2. Mean streamwise base velocity, U{y), obtained in 2D simulations of the uncon- 
trolled Poiseuille flow with = 63.25, x, and controlled flows subject to: (a) DTWs with 
□ , (c = 5, cj^ = 2, a = 0.035); o, (c = 5, cj^ = 2, a = 0.05); 0, (c = 5, cj^ = 2, a = 0.125); 
and (b) UTWs with <, (c = -2, uo^ = 0.5, a = 0.015); V, (c = -2, cj^. = 0.5, a = 0.05); 
A, (c = -2, cj^ = 0.5, a = 0.125). 



As shown by Hoepffner fc Fukagata| ( |2009 ), compared to the uncontrolled laminar flow, 
the nominal flux is reduced (increased) by DTWs (UTWs); according to (2.3), this results 
in larger (smaller) nominal drag coefficients, respectively. 

The above results suggest that properly chosen traveling waves can exhibit increased 
flux compared to the uncontrolled flow. For flxed pressure gradient, this results in pro- 
duction of a driving power 

where Ub,c and Ub,u denote the flux of the controlled and uncontrolled flows. The nor- 
malized produced power %Ilprod is expressed as a percentage of the power spent to drive 
the uncontrolled flow, = —P^ Ub,u i'^L^Lz)^ 

%Ilprod = 100(/7b,c - Ub,u)/Ub,u- 

On the other hand, the input power required for maintaining the traveling waves is 
obtained from ( Currie||2003 ) 

and the normalized required power %Ilreq is expressed as 

VP\ . - VP\ _ 

%Ureq = 100 — 



-2P,Ub,u 

In order to assess the efficacy of traveling waves for controlling transitional flows, the 
control net power is deflned as the difference between the produced and required pow- 
ers dQuadrio fc Ricco||2004D 

where %Ilnet signifles how much net power is gained (positive %Ilnet) or lost (negative 
%Ilnet) in the controlled flow as a percentage of the power spent to drive the uncontrolled 
flow. 
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Case 


c 




a 


Ub,n 


10' Cf,N 




%Iireq 













0.6667 


4.5000 











1 


5 


2 


0.035 


0.6428 


4.8404 


-3.58 


16.64 


-20.22 


2 


5 


2 


0.050 


0.6215 


5.1778 


-6.77 


31.74 


-38.51 


3 


5 


2 


0.125 


0.4821 


8.6050 


-27.69 


136.50 


-164.19 


4 


-2 


0.5 


0.015 


0.6703 


4.4513 


2.70 


5.46 


-2.76 


5 


-2 


0.5 


0.050 


0.7791 


3.2949 


16.86 


37.69 


-20.83 


6 


-2 


0.5 


0.125 


1.0133 


1.9478 


51.99 


145.05 


-93.06 



Table 2. Nominal results in Poiseuille flow with Rr — 63.25. The nominal flux, U b,n , and 
skin-friction drag coefficient, C/,Ar, are computed using the base flow described in § |2.3| The 
produced power, YoUprod, required power, %Ilreq, and net power, %,Ilnet, are normalized by the 
power required to drive the uncontrolled flow. The produced and net powers are computed with 
respect to the laminar uncontrolled flow. 



The nominal efficiency of the selected streamwise traveling waves in 2D flows, i.e. in 
the absence of velocity fluctuations, is shown in table |2j Note that the nominal net power 
is negative for all controlled 2D simulations. This is in agreement with a recent study 
of Hoepffner & Fukagata ( 2009 ) where it was shown that the net power required to drive 
a flow with wall transpiration is always larger than in the standard pressure gradient 
type of actuation. 



3. Avoidance/promotion of turbulence by streamwise traveling waves 

In Part 1 it was shown that a positive net efficiency can be achieved in a situation where 
the controlled flow stays laminar but the uncontrolled flow becomes turbulent. Whether 
the controlled flow can remain laminar depends on velocity fluctuations around the mod- 
ified base flow. In this section, we study the influence of streamwise traveling waves on 
the dynamics and the control net efficiency. This problem is addressed by simulating a 3D 
channel flow with initial perturbations which are superimposed on the base velocity in- 
duced by the wall actuation. Depending on the kinetic energy of the initial condition, we 
distinguish three cases: (i) both the uncontrolled and properly designed controlled flows 
remain laminar {small initial energy); (ii) the uncontrolled flow becomes turbulent, while 
the controlled flow stays laminar for the appropriate choice of traveling wave parameters 
{moderate initial energy); and (iii) both the uncontrolled and controlled flows become 
turbulent for selected traveling wave parameters {large initial energy). Our simulations 
indicate, however, that poorly designed traveling waves can promote turbulence even for 
initial conditions for which the uncontrolled flow stays laminar. It was demonstrated in 
Part 1 that properly designed DTWs are capable of significantly reducing receptivity of 
velocity fluctuations which makes them well-suited for preventing transition; on the other 
hand, compared to the uncontrolled flow, the velocity fluctuations around the UTWs at 
best exhibit similar receptivity to background disturbances. Following Part 1, we present 
our main results for DTWs with (c = 5, co^^^ = 2); these results are compared to UTWs 
with (c = —2, ujx = 0.5) (as selected in Min et al. ( |2QQ6| )). In both cases, three wave 
amplitudes are selected (cf. table [T]). 

The 3D simulations, which are summarized in table [3j confirm and complement the 
theoretical predictions of Part 1 at two levels. At the level of controlling the onset of 



turbulence, we illustrate in § 3.1 that the UTWs increase receptivity of velocity fluctu- 



ations and promote turbulence even for initial perturbations for which the uncontrolled 
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Initial Energy Case c Ux a 10^ Cf Vollprod %^req %^net 



Small 












4.5002 











2 


5 


2 


0.050 


5.1778 


-6.77 


31.77 


-38.54 


4 


-2 


0.5 


0.015 


4.3204 


-1.54 


5.14 


-3.60 


5 


-2 


0.5 


0.050 


5.9426 


-16.52 


23.22 


-39.74 


6 


-2 


0.5 


0.125 


3.6853 


12.20 


108.41 


-96.21 



Moderate 












10.3000 











1 


5 


2 


0.035 


4.9244 


52.07 


26.44 


25.63 


2 


5 


2 


0.050 


5.2273 


47.35 


50.40 


-3.05 


4 


-2 


0.5 


0.015 


8.7866 


11.36 


4.53 


6.83 


5 


-2 


0.5 


0.050 


6.7406 


31.15 


41.96 


-10.81 


6 


-2 


0.5 


0.125 


3.9264 


77.03 


155.80 


-78.77 



Large 












11.2000 











2 


5 


2 


0.050 


11.9000 


-3.37 


47.90 


-51.27 


3 


5 


2 


0.125 


12.1000 


-11.31 


196.89 


-208.20 


5 


-2 


0.5 


0.050 


7.4438 


13.68 


34.19 


-20.51 


6 


-2 


0.5 


0.125 


3.9872 


57.75 


142.92 


-85.17 



Table 3. Results of 3D simulations in Poiseuille flow with Rr = 63.25 for initial conditions of 
small, moderate, and large energy (respectively, ^(0) = 2.25 x 10~^, E{0) — 5.0625 x 10~^, 
and ^(0) = 2.5 x 10~^). The values of C/, Voliprod, and Xllnet correspond to t = 1000. 

For small initial energy, the produced and net powers are computed with respect to laminar 
uncontrolled flow; for moderate and large initial energies, they are computed with respect to 
turbulent uncontrolled flow. 



flow stays laminar. In contrast, the DTWs can prevent transition even in the presence of 



initial conditions with moderate and large energy (cf. § |3.2| and § |3.3[ ) . At the level of net 
power efficiency, it is first shown in § |3.1| that the net power is negative when the uncon- 
trolled flow stays laminar. However, for the uncontrolled flow that becomes turbulent, we 



demonstrate that the DTWs can result in a positive net efficiency. As discussed in § 3.2 
and § |3.3[ the positive net efficiency is achieved if the required power for maintaining 
the laminar DTW is less than the produced power. In addition, in § |3.3[ we highlight 
an important trade-off that limits the advantages of DTWs in controlling the onset of 
turbulence in flows subject to large initial conditions. Namely, we show that in this case 
preventing transition by DTWs requires a large input power that results in a negative 



efficiency. Our simulations in § |3. 2 [ reveal that although UTWs become turbulent, a pos- 
itive net efficiency can be achieved for small enough wave amplitudes. For the initial 
conditions with moderate energy, we further point out that the achievable positive net 
efficiency for UTWs is much smaller than for the DTWs that sustain the laminar flow 
(cf. §|3l). 



3.1. Small initial energy 

We first consider the initial perturbations with small kinetic energy, £^(0) = 2.25 x 10~^, 
which cannot trigger turbulence in flow with no control. Our simulations show that the 
DTWs selected in Part 1 of this study improve transient response of the velocity fluctu- 



ations; on the contrary, the UTWs considered in Min et al. (2006) lead to deterioration 
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Figure 3. Energy of the velocity fluctuations, E{t), for the initial condition with small energy: 
(a) X, uncontrolled; o, a DTW with (c = 5, cj^. = 2, a = 0.05); and (b) UTWs with <, (c = -2, 
= 0.5, a = 0.015); V, (c = -2, Ua: = 0.5, a = 0.05); A, (c = -2, cj,, = 0.5, a = 0.125). 



of the transient response and, consequently, promote turbulence. Since the uncontrolled 
flow stays laminar, both DTWs and UTWs lead to the negative net efficiency. 
The energy of velocity fluctuations is given by 



where ft = 2LxLz is the volume of the computational box. Figure |3] shows the fluctua- 
tions' kinetic energy as a function of time for the uncontrolled flow and controlled flows 
subject to a DTW with (c = 5, cj^ = 2, a = 0.05) and three UTWs with (c = -2, 
Ux = 0.5, a = {0.015, 0.05, 0.125}). As evident from figure 3(a) , the energy of the uncon- 
trolled flow exhibits a transient growth followed by an exponential decay to zero (i.e., to 
the laminar flow). We see that a DTW moves the transient response peak to a smaller 
time, which is about half the time at which peak of E{t) in the uncontrolled flow takes 
place. Furthermore, maximal transient growth of the uncontrolled flow is reduced by 
approximately 2.5 times, and a much faster disappearance of the velocity fluctuations is 



achieved. On the other hand, figure 3(b) clearly exhibits the negative influence of the 
UTWs on a transient response. In particular, the two UTWs with larger amplitudes 
significantly increase the energy of velocity fluctuations. We note that the fluctuations' 
kinetic energy in a flow subject to a UTW with an amplitude as small as a = 0.015 at 
t = 1000 is already about two orders of magnitude larger than the maximal transient 
growth of the flow with no control. 



Figure 4(a) shows the skin- friction drag coefficient. 



Cf{t) 



2r 
^ B 



1 



Rc U% 



dU du 
dy dy 



dU du\ 
dy dy) 



2/ = ! 



as a function of time for the traveling waves considered in figure [3] Here, denotes the 
non-dimensional average wall-shear stress and 



{U{y) + u{y,t)) dy, 



is the total bulk flux. Since both the uncontrolled flow and the flow subject to a DTW 
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Figure 4. (a) Skin- friction drag coefficient, C/; (b) normalized required power, %nreg; (c) nor- 
malized produced power, %Ilprod'i and (d) normalized net power, %Ilnet, for the initial condition 
with small energy: x, uncontrolled flow; o, DTW with (c = 5, cj^c = 2, a = 0.05); and UTWs 
with <, (c -2, 0.5, a 0.015); V, (c -2, cj^ 0.5, a 0.05); A, (c = -2, cj^ 0.5, 

a = 0.125). 



stay laminar, their steady-state drag coefficients agree with the nominal values computed 
in the absence of velocity fluctuations (cf. tables |2] and [3| . On the other hand, the drag 
coefficients of the UTWs that become turbulent are about twice the values predicted 
using the base flow analysis. The large amplification of velocity fluctuations by UTWs is 
responsible for this increase. The velocity fluctuations in the UTW with a = 0.015 are 
not amplified enough to have a pronounced effect on the drag coefficient. Furthermore, 
the drag coefficients for t he UTWs with (c = -2, uj^ = 0.5, a = {0.05, 0.125}) at t = 1000 
agree with the results of'Min e t al.\ ( 2QQ6[ ) computed for the fully developed turbulent 
channel flow. This indicates that the UTWs with larger amplitudes in our simulations 
have transitioned to turbulence. The above results confirm the theoretical prediction of 
Part 1 where it is shown that the UTWs are poor candidates for controlling the onset of 
turbulence for they increase receptivity relative to the uncontrolled flow. 

The normalized required, produced, and net powers for the initial conditions with small 



kinetic energy are shown in figures 4(b) - 4(d) Note that the normalized net power for all 



traveling waves is negative (cf. figure 4(d)[ ). This confirms the prediction of Part 1 that 
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the net power is negative whenever the uncontrohed flow stays laminar. It is noteworthy 
that the UTW with a = 0.125 has a negative net power despite its significantly smaller 
drag coefficient compared to the laminar uncontrolled flow. As evident from figures |4(b)| 
and |4(c)[ this is because the required power for maintaining this UTW is much larger 
than the power produced by reducing drag. The above results agree with the studies 
of |Bewley ( 2009 ) and Fukagata et al ( 2009 ) where it was established that the net cost to 



drive a flow by any transpiration-based strategy is larger than in the uncontrolled laminar 
flow. Therefore, aiming for sub-laminar drag may not be advantageous from efficiency 
point of view. Instead, one can design control strategies that yield smaller drag than the 
uncontrolled turbulent flow and provide positive net power balance (cf. § |3.2[). 



3.2. Moderate initial energy 

We next consider the velocity fluctuations with moderate initial energy, E{0) = 5.0625 x 
10~^. This selection illustrates a situation where the initial conditions are large enough 
to trigger turbulence in the uncontrolled flow but small enough to allow the properly 
chosen DTWs to maintain the laminar flow and achieve positive net power balance. As 
shown in § |3.1[ the UTWs trigger turbulence even for the initial conditions whose kinetic 
energy is about 200 times smaller than the value considered here. 

The energy of the velocity fluctuations and drag coefficients as a function of time for the 
uncontrolled flow and UTWs are shown in flgures 5(a) and |5(b)] Figure 5(a)| indicates 



that the kinetic energy of the uncontrolled flow and the flow subject to UTWs with 
(c = —2, (jOx = 0.5, a = {0.015,0.05,0.125}) is increased by orders of magnitude which 
eventually results in transition to turbulence. This large energy ampliflcation of UTWs is 
captured by the linear analysis around the laminar base flows in Part 1. As evident from 
flgure |5(b)[ the large fluctuations' energy in both the uncontrolled flow and in UTWs 
yields much larger drag coefficients compared to the nominal values reported in table |2j 
In addition, figure [5(b)] is in agreement with |Min et al. (2006|) where it was shown that 



the skin-friction drag coefficients of the UTWs are smaller than in the uncontrolled flow 
that becomes turbulent, and that the UTW with (c = —2, uj^ = 0.5, a = 0.125) achieves 
a sub-laminar drag. The normalized required and net powers for the UTWs are shown 
in flgures 5(c) and |5(d)[ Note that the required power for maintaining the UTW with 



a = 0.125 (which yields sub-laminar drag) is so large that it results in a negative net 
power balance (cf. flgure [S^d) ). On the other hand, the UTW with a = 0.015 is capable 



of producing a small positive net power for two main reasons: (i) it has a smaller drag 
coefficient than the uncontrolled turbulent fiow (although it becomes turbulent itself); 
and (ii) it requires a much smaller power compared to the UTW with a = 0.125. 

The fluctuations' kinetic energy and skin-friction drag coefficient for the DTWs are 
shown in figures |6(a)| and |6(b)[ Figure 6(a) shows that the DTWs with (c = 5, 



2, a = {0.035,0.05}) significantly weaken intensity of the velocity fiuctuations, thereby 
facilitating maintenance of the laminar fiow. From figure [6(b)] we also see that the small 
transient growth of fiuctuations' kinetic energy results in a small transient increase in 
the drag coefficients which eventually decay to their nominal values reported in table [2] 
Even though these drag coefficients are larger than in the uncontrolled laminar flow, they 
are still approximately two times smaller than in the uncontrolled flow that becomes 
turbulent (cf. table [3|. 

The normalized produced, required, and net powers for DTWs are shown in flgures 6(c)| 



and |6(d)[ As can be seen from flgure [6(c)| the normalized produced power for the DTWs 
is positive by virtue of the fact that the uncontrolled flow becomes turbulent while 
the controlled flows stay laminar. Figure |6(d)| shows that the DTW with a = 0.035 
(respectively, a — 0.05) has a positive (respectively, negative) net power balance. The 
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Figure 5. (a) Energy of the velocity fluctuations, E(t); (b) skin- friction drag coefficient, Cf{t); 
(c) normalized required power, %Ilreq; and (d) normalized net power, %Ilnet, for the initial con- 
dition with moderate energy: x, uncontrolled; and UTWs with <, (c = —2,Ux — 0.5, a — 0.015); 
V, (c= -2,cj^ ^ 0.5, a = 0.05); A, (c = -2,cj^ = 0.5, a = 0.125). 



reason for this is twofold: first, as evident from figure |6(c)[ the DTW with lar ger a results 
in a smaller produced power since it induces a larger negative nominal bulk flux than 
the DTW with smaller a; and second, the required power to maintain the DTW with 
larger a is bigger than in the DTW with smaller a. Furthermore, at t = 1000, the DTW 
with a = 0.035 has a larger net power than the UTW with a = 0.015 {%Ilnet = 25.63 
vs. %Ilnet = 6.83; cf. table [3|. This is because the DTW with a = 0.035, in contrast to 
the UTW with a = 0.015, remains laminar and produces a much larger power than it 
requires. 

In summary, the results of this section highlight an important trade-off that needs to 
be taken into account when designing the traveling waves. Large amplitudes of properly 
designed downstream waves yield larger receptivity reduction which is desirable for con- 
trolling the onset of turbulence. However, this is accompanied by an increase in drag 
coefficient and required control power. Thus, to maximize net efficiency, it is advanta- 
geous to select the smallest possible amplitude of wall-actuation that can maintain the 
laminar flow. 
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Figure 6. (a) Energy of the velocity fluctuations, E{t)] (b) skin- friction drag coefficient, Cf{t); 
(c) normalized required power, %Ilreq (solid), normalized produced power, %Ilprod (dashed); 
and (d) normalized net power, %nnet, for the initial condition with moderate energy: DTWs 
with □, (c = 5,cJx = 2,a = 0.035); o, (c = 5, cj^ = 2, a = 0.05). 



3.3. Large initial energy 



Section |3.2| illustrates capability of properly designed DTWs to maintain the laminar 
flow in the presence of initial conditions that induce transition in the uncontrolled flow. 
In this section, we demonstrate that, as the energy of the initial perturbation increases, 
a DTW with larger amplitude is needed to prevent transition. Our results confirm the 
prediction made in Part 1 that maintaining a laminar flow with a larger DTW amplitude 
comes at the expense of introducing a negative net power balance. 

Simulations in this section are done for the initial condition with large kinetic energy, 
E{0) = 2.5 X 10~^. The time evolution of the fluctuations' energy for a pair of DTWs with 
(c = 5, cja^ = 2, a = {0.05, 0.125}) is shown in figure [7(a)] The uncontrolled flow becomes 
turbulent and exhibits similar trends in the evolution of E{t) as the corresponding flow 
initiated with moderate energy perturbations (cf. figures [5(a)] and [7(a) [ ). On the other 
hand, figure [7(a)] shows that the DTW with a = 0.05 is not capable of maintaining the 
laminar flow; in comparison, the same set of control parameters prevented transition for 
the perturbations of moderate initial energy (cf. figures 6(a) and |7(a)] ). Conversely, the 




Figure 7. (a) Energy of the velocity fluctuations, E(t); (b) skin- friction drag coefficient, Cf{t); 

(c) normalized required power, %Ilreq (solid), normalized produced power, %Ilprod (dashed); and 

(d) normalized net power, %,Ilnet, for the initial condition with large energy: x, uncontrolled; 
DTWs with o, (c= 5,cJx = 2,a = 0.05); and 0, (c = 5, cj^ = 2, a = 0.125). 



DTW with a = 0.125 remains laminar even though E(t) transiently reaches about half 
the energy of the turbulent uncontrolled flow. Therefore, the DTWs with frequency and 
speed selected in Part 1 and sufficiently large amplitudes are capable of maintaining the 
laminar flow even in the presence of large initial perturbations. 

Figure |7(b)| shows the skin-friction drag coefficients for the flows considered in fig- 



ure 



7(a)[ For the DTW with a = 0.05 the steady-state value of Cf is given by C/ 



11.9 X 10-^, which is a slightly larger value than in the turbulent uncontrolled flow, 
Cf = 11.2 X 10"^ (cf. table [3|. We note that the drag coefficient of the DTW that stays 
laminar initially reaches values that are about 50 % larger than in the uncontrolled flow; 
after this initial increase, Cf{t) then gradually decays to the va lue predicted using the 
base flow analysis, C/ = 8.6 x 10~^ (cf. table [2|. Figures [7(^ and [7(d)] show the nor- 
malized required, produced, and net powers for the initial condition with large kinetic 
energy. As evident from figure |7(d)[ the net power balance is negative for all considered 
flows. The DTW with a = 0.05 becomes turbulent, and it has a larger drag coefficient 
than the uncontrolled flow which consequently leads to negative produced and net pow- 
ers. Moreover, even though the DTW with a = 0.125 can sustain laminar flow, its net 
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Figure 8. (a) and (b) Production, PEif) (solid), and dissipation, DEit) (dashed), of kinetic 
energy in Poiseuille flow with Rr — 63.25 for the initial condition with small energy: (a) x, 
uncontrolled; o, DTW with (c = 5, cj^: = 2, a = 0.05); and (b) V, UTW with (c = -2, cj^ = 0.5, 
a = 0.05); A, UTW with a = 0.125. (c) PE{t) + DE{t) for the UTW with (c = -2, cj^ = 0.5, 
a = 0.05). 



power balance is very poor. There are two main reasons for the lack of efficiency of this 
control strategy: first, its nominal drag coefficient is significantly larger than in a DTW 
with smaller amplitudes which consequently yields very small produced power (at larger 
times not shown in figure 7(c)); and second, a prohibitively large power is required to 
maintain this large amplitude DTW. 

The results of this section show that preventing transition by DTWs in the presence 
of large initial conditions comes at the expense of large negative net power balance. We 
also highlight that in the presence of large initial perturbations (or, equivalent ly, at large 
Reynolds numbers), transition to turbulence may be inevitable. Furthermore, the results 
of § |3.2| show that the UTWs may reduce the turbulent skin-friction drag and achieve 
positive net efficiency. The approach used in Part 1 considers dynamics of fluctuations 
around laminar flows and, thus, it cannot be used for explaining the positive efficiency 
of the UTWs that become turbulent. 
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3.4. Energy amplification mechanisms 

The Reynolds-Orr equation can be used to quantify the evolution of fluctuations' kinetic 
energy around a given base flow ( Schmid fc Henningson|2QQT ). In this section, we use this 



equation to examine mechanisms that contribute to production and dissipation of kinetic 
energy in flows subject to traveling waves. For base velocity, U5 = {U (x, t), V{x^ y^ t)^0)^ 
the evolution of the energy of velocity fluctuations, E{t), is determined by 

PE{t) = -^ij i^^^v + ^^^^ + ^^^v + u^U^)<i^, (3.1) 

Here, Pe represents the production of kinetic energy and is associated with the work of 
the Reynolds stresses on the base shear, whereas De accounts for viscous dissipation. 

We confine our attention to the simulations for initial conditions with small energy. 
This situation is convenient for explaining why the DTWs exhibit improved transient 
behavior compared to the laminar uncontrolled flow while the UTWs promote transition 



to turbulence. Figure 8(a) shows production and dissipation terms for the uncontrolled 
flow and for the flow subject to a DTW with (c = 5, Ux = 2, a = 0.05). For the 
uncontrolled flow, Pe is always positive, De is always negative, and they both decay 
to zero at large times. On the contrary, the production term for the DTW becomes 
negative for SO < t < 220. We see that, at early times, Pe and De for the DTW follow 
their uncontrolled flow counterparts. However, after this initial period, they decay more 
rapidly to zero. These results confirm the prediction of Part 1 that the DTWs reduce 
the production of kinetic energy. In contrast, figure |8(b)| shows that the UTWs with 
(c = —2, Ux = 0.5, a = {0.05,0.125}) increase both Pe and De by about four orders of 
magnitude compared to the values reported in figure |8(a)[ This verifies the theoretical 
prediction of Part 1 that the UTWs increase the production of kinetic energy. Moreover, 
production dominates dissipation transiently, thereby inducing large growth of kinetic 
energy observed in figure |3(b)[ For the UTW with a = 0.05, this is further illustrated 



in figure 8(c) by showing that Pe accumulates more energy than De dissipates (i.e., 
the area under the curve in figure 8(c)| is positive). We also note that, in the above 



simulations, the work of Reynolds stress uv on the base shear Uy dominates the other 
energy production terms. Furthermore, our results show that the wall-normal diffusion 
of u is responsible for the largest viscous energy dissipation. 



3.5. Flow visualization 

In sections |3.1| - 13.3[ transition was identified by examining fluctuations' kinetic energy 
and skin-friction drag coefficients. Large levels of sustained kinetic energy and substantial 
increase in drag coefficients (compared to base flows) were used as indicators of transition. 
Here, we use flow visualization to identify coherent structures in both the uncontrolled 
and controlled flows. 

The onset of turbulence in a bypass transition is usually characterized by formation 
of streamwise streaks and their subsequent break-down. For the initial condition with 
moderate energy, figure |9] shows the streamwise velocity fluctuations at = —0.5557 
(7/+ = 28.11 in wall units) for the uncontrolled flow and flows subject to a UTW with 
(c -2, Ux 0.5, a 0.05) and a DTW with (c = 5, cj^ = 2, a = 0.05). Cle arly, 
the initial perturbations evolve into streamwise streaks in all three flows (cf. figures [9 (a) | 
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Figure 9. Streamwise velocity fluctuations, u{x,z), at y = —0.5557 {y~^ = 28.11), (a) - (c) 
t = 50, and (d) - (f) t = 120 for initial condition with moderate energy: uncontrolled flow; UTW 
with (c -2, cj^ 0.5, a 0.05); and DTW with (c = 5, cj^ = 2, a = 0.05). 



9(c)| for t = 50). At t = 120, the growth of velocity fluctuations results in a break-down 



of the streaks both in the flow with no control and in the flow subject to UTWs (cf. 
figures [9(d)| and 9(e)). We see that the streaks evolve into complex flow patterns much 
faster in the latter case. For the UTWs, the streak distortion occurs as early as t = 50 and 
a broad range of spatial scales is observed at t = 120. On the contrary, figure [9(f)] shows 
that, at t = 120, the DTWs have reduced the magnitude of velocity fluctuations to about 
half the value at t = 50, thereby weakening intensity of the streaks and maintaining the 
laminar flow. 



4. Relaminarization by downstream waves 

Thus far we have shown that properly designed DTWs represent an effective means 
for controlling the onset of turbulence. In this section, we demonstrate that the DTWs 
designed in Part 1 can also relaminarize fully developed turbulent flows. Since the lifetime 
of turbulence depends on the Reynolds number (Brosa |1989[|Grossmann|2QQQ||Hof et al 



2QQ6| [Borrero-Echeverry et fl/.pQlO] ), we examine turbulent flows with Rc = 2000 (i.e.. 



Rr ^ 63.25) and Rc = 4300 (i.e., Rr ^ 92.80). 

The numerical scheme described in § |2.2| is used to simulate the turbulent flows. For 
Rc = 4300, the number of grid points in the streamwise, wall-normal, and spanwise 
directions is increased to 80 x 97x 80. The velocity field is initialized with a fully developed 
turbulent flow obtained in the absence of control. The surface blowing and suction that 
generates DTWs is then introduced (at t = 0), and the kinetic energy and drag coefficients 
are computed at each time step. 

The fluctuations' kinetic energy for the uncontrolled flow and for the flows subject 
to DTWs with (c = 5, cj^ = 2, a = {0.035, 0.05, 0.125}) at Rc = 2000 are shown in 
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Figure 10. Energy of velocity fluctuations around base flows of § |2.1| Simulations are ini- 
tiated by a fully developed turbulent flow with Rc = 2000: x, uncontrolled; DTWs with 
□ , (c = 5,cJa. = 2,a = 0.035); o, (c = 5, cj^. = 2, a = 0.05); and 0, (c = 5, cJo. = 2, a = 0.125). 
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Figure 11. (a) Energy of velocity fluctuations around base flows of § 2.1 E{t); and (b) 
skin- friction drag coefficient, Cf{t). Simulations are initiated by a fully developed turbu- 
lent flow with Rc = 4300: x, uncontrolled; DTWs with □, (c = 5,cJa: = 2,a = 0.035); 
o, (c = 5,cj^ = 2,a = 0.05); and 0, (c = 5, cj^ = 2, a = 0.125). 



figure 10 The energy of velocity fluctuations around base flows U5 of § |2.1| (parabola for 
flow with no control; traveling waves for flow with control) are shown; relaminarization 
occurs when the energy of velocity fluctuations converges to zero. Clearly, large levels 
of fluctuations in the flow with no control are maintained up until t ~ 1100. After this 
time instant, however, velocity fluctuations exhibit gradual decay. On the other hand, 
fluctuations in flows subject to DTWs start decaying much earlier, thereby indicating 
that the lifetime of turbulence is reduced by surface blowing and suction. Relative to the 
uncontrolled flow, the fluctuations' kinetic energy for the DTWs considered here converge 
much faster to zero. We also see that the rate of decay increases as the wave amplitude 
gets larger. 

We next consider a turbulent flow with Rc = 4300. The energy of velocity fluctua- 
tions around base flows U5 of § 2.1 is shown in figure pT(a)[ In both the uncontrolled 
flow and the flows subject to the DTWs with (c = 5, a;^^ = 2, a = {0.035,0.05}) the 
energy oscillates around large values that identify turbulent flow. This indicates that the 
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Figure 12. Mean velocity, U{y), and streamwise velocity (color plots) ai y = —0.7518. Simula- 
tions are initiated by a fully developed turbulent flow with Rc = 4300: x,t = 0;«, t=100;*, 
t = 900. The dashed line identifies the laminar mean velocity induced by the DTW with (c = 5, 
cj^ = 2, a = 0.125). 



DTWs with smaller amplitudes cannot eliminate turbulence. On the contrary, the DTW 
with a = 0.125 reduces the energy of velocity fluctuations, thereby relaminarizing the 



flow. Figure |ll(b) shows that the skin-friction drag coefficient for the uncontrolled flow 



and for the DTWs with smaller amplitudes is approximately constant throughout the 
simulation. On the other hand, owing to relaminarization, the drag coefficient for the 
DTW with a = 0.125 is smaller than that of the uncontrolled turbulent flow. However, 
relaminarization comes at the expense of poor net efficiency. This is because of the large 
required power (i.e., high cost of control), which reduces the appeal of using DTWs for 
control of turbulent flows. 

Figure 12 shows the mean velocity, U{y)^ at three time instants in the flow subject 
to the DTW with (c = 5, cj^^ = 2, a = 0.125). The instantaneous values of stream- 
wise velocity in the (x, z)-plane at = —0.7518 are also shown. As time advances, the 
initial turbulent mean velocity at Rc = 4300 moves towards the laminar mean velocity 
induced by the surface blowing and suction (dashed line). The color plots illustrate how 
the initial turbulent flow evolves into the DTW laminar base flow. We conjecture that 
receptivity reduction is important not only for controlling the onset of turbulence but 
also for relaminarization of fully developed flows. As outlined in Part 1, explaining the 
effect of traveling waves on turbulent flows requires additional control-oriented modeling 
and further scrutiny. 



5. Concluding remarks 

This work, along with a companion paper (Moarref & Jovanovic 2010), represents a 
continuation of recent efforts (Jovanovic, Moarref & You 2006 Jovanovic|2008 ) to develop 
a model-based paradigm for design of sensor less flow control strategies in wall-bounded 
shear flows. Direct numerical simulations are used to complement and verify theoretical 
predictions of Moarref fc Jovanovic | ( [2010 ) , where receptivity analysis of the linearized NS 
equations was used to design small amplitude traveling waves. We have shown that per- 
turbation analysis (in the wave amplitude) represents a powerful simulation-free method 
for predicting full-scale phenomena and controlling the onset of turbulence. 
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Simulations of nonlinear flow dynamics have demonstrated that the DTWs, designed 
in Part 1, can maintain laminar flow and achieve positive net efficiency. In contrast, the 
UTWs promote turbulence even with the initial conditions for which the uncontrolled 
flow remains laminar. Our analysis of the Reynolds-Orr equation shows that, compared 
to the uncontrolled flow, the DTWs (UTWs) reduce (increase) the production of kinetic 
energy. 



We have also examined the effects of DTWs on fully developed turbulent flows at low 
Reynolds numbers. It turns out that the DTWs with speed and frequency selected in 
Part 1 and large enough amplitudes can eliminate turbulence (i.e., relaminarize the flow). 
We also note that, in spite of promoting turbulence, the UTWs may still achieve smaller 
drag coefficients compared to the uncontrolled flow. By increasing the UTW amplitude, 
even sub-laminar drag can be attained ( Min et a/.|2QQ6 ). It is to be noted, however, that 
large wave amplitudes introduce poor net efficiency in flows subject to a fixed pressure 
gradient. Nevertheless, these traveling waves may still be utilized when the primal interest 
is to eliminate turbulence (with DTWs) or reduce the skin-friction drag (with UTWs) 
irrespective of the cost of control. 



All simulations in the present study are enforced by a fixed pressure gradient, as op- 



posed to the constant mass flux simulations of Min et al. (2006). This is consistent with 



Part 1 where receptivity analysis was done for flows driven by a fixed pressure gradient. 
Even though these setups are equivalent in steady flows (|Scotti & Piomelli 2001), they 



can exhibit fundamentally different behavior in unsteady flows. For example, the two 
simulations may possess different stability characteristics and yield structurally different 
solutions for near- wall turbulence ( Waleffe|2001 Jimenez et a/.|20Q5 ). Moreover, Jimenez 
et al\ ([2005) remarked that the dynamical properties of the two simulations can sig- 
niflcantly differ in small computational domains. Also, Kerswell ( |2005J suggested that 
specifying the simulation type comes next to deflning the boundary conditions. 



In order to examine the effect of simulation type on transition, skin-friction drag coef- 
flcient, and control net efficiency, we have repeated some of the simulations by adjusting 
the pressure gradient to maintain a constant mass fiux. Our results reveal that regard- 
less of the simulation type, the DTWs designed in Part 1 are effective in preventing 
transition while the UTWs promote turbulence. Moreover, the steady-state skin-friction 
drag coefficients are almost identical in both cases. However, the control net efficiency 
depends significantly on the simulation type. This is because of the difference in the def- 
inition of the produced power: in the fixed pressure gradient setup, the produced power 
is captured by the difference between the bulk fluxes in the uncontrolled and controlled 
flows; in the constant mass flux setup, the produced power is determined by the differ- 
ence between the driving pressure gradients in the uncontrolled and controlled flows. It 
turns out that the produced power is larger in the constant mass flux simulation than in 
the flxed pressure gradient simulation, whereas the required power remains almost un- 
changed. Consequently, both DTWs and UTWs have larger efficiency in constant mass 
flux simulations. For example, in fixed pressure gradient setup of the present study, the 
UTWs with (c = —2^uOx = 0.05 a = {0.05, 0.125}) have negative efficiency. The efficiency 
of these UTWs is positive, however, in constant mass fiux simulations ( Min et ar||2006 ). 
Our ongoing effort is directed towards understanding the reason behind this disagree- 
ment which may be ultimately related to the fundamental difference between these two 
types of simulations. 
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